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We study the recurrence coefficients of the orthogonal polynomials with respect 
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i> ■ 1 Introduction 

o 

■ 1.1 Orthogonal polynomials 

^ ! It is well-known (see, for instance, |3l |T^) that a sequence of orthonormal polynomials 

(Pn)ngN, i-e.. 



Pra{x)Pn{x) dfi{x) = Sm,n, (1-1) 

where is a positive measure with support on the real line and 6m,n is the Kronecker 
delta, satisfies a three-term recurrence relation 

xpn{x) = a„+ip„+i(x) + hnPn{x) + a„p„_i(a;) (1.2) 

with the recurrence coefficients given by the following integrals 

xpn{x)pn^i{x)diJ,{x), 6„ = xpl{x) dfi{x) . (1.3) 



Here it is assumed that p_i = 0. 

We will consider discrete orthogonal polynomials for which the weight is supported on 
an equidistant lattice {hn -|- no| n G A C Z} with parameters h and uq. In case of the 
lattice No = {0, 1, . . .} the orthogonality condition (11. ip reads 

oo 

^Pm{k)pn{k) W{k) = 6m,n- (1-4) 
k=0 



The examples of (classical) discrete orthogonal polynomials on an equidistant lattice in- 
clude the Meixner and Charlier polynomials with the lattice Nq and the Krawtchouk and 
Hahn polynomials (A^ + 1 polynomials orthogonal on {0, 1, . . . , N}). For more infor- 
mation about discrete (classical) orthogonal polynomials, we refer to [191 ES] (see also 
[21 ini [in] for the definitions of classical weights and their semi-classical extensions and 
further references). 

A useful characterization of classical polynomials is the Pearson equation 

[cr{x)w{x)]' = t{x)w{x), 

where a and r are polynomials satisfying degcr < 2 and degr = 1. In case of discrete 
polynomials on a linear lattice the Pearson equation reads 

V[a{x)w{x)] = t{x)w{x), 

where V is the backward difference operator 

V/(x) = /(x)-/(x-l). 

Semi-classical orthogonal polynomials are defined as orthogonal polynomials for which the 
weight satisfies a Pearson equation for which degcr > 2 or degr 7^ 1 (see [HI EH]). The 
recurrence coefficients in the three-term recurrence relation (11. 2p for classical orthogonal 
polynomials (e.g., Hermite, Laguerre, Jacobi, Charlier, Meixner, Krawtchouk and others) 
can be found explicitly in contrast to non-classical weights. The recurrence coefficients of 
semi-classical weights usually obey nonlinear recurrence relations, which, in many cases, 
can be identified as discrete Painleve equations; see for instance [31 [22] and the references 
therein. Moreover, when viewed as functions of one of the parameters in the weight, the 
recurrence coefficients satisfy (continuous) Painleve equations (see [21 [3 [10] in case of the 
discrete semi-classical Meixner and Charlier weights and fLl] and the references therein 
for other semi-classical continuous weights). This paper gives another example of the 
connection between discrete orthogonal polynomials and the Painleve equations. 

1.2 Statement of the results 

The classical Krawtchouk polynomials for each N and n < N are defined by 

Kn{x) = 2F,(^-n,-x--N;^y pe(0,l), (1.5) 

where 

.F.(..6;e;.) = f Mii^.-, 
with the Pochhammer symbol (a)„ being defined by 

n-l 

(a)„ = ]^(a + j) = a{a + 1) ■ ■ ■ (a + - 1), 

j=0 

is the hypergeometric function; cf. [H Chapter 15]. They satisfy orthogonality conditions 
with respect to the binomial weight 

w{k) = (^]p\i - pf-', ke{0,i,..., N}. 
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The polynomials satisfy fll.2p for n = 0, . . . , N — 1 with K_i = and = 0. The 
Krawtchouk polynomials are, in fact, the Meixner polynomials 



with P = —N and i = p/{p ~ 1). The recurrence coefficients for the orthonormal 
Krawtchouk polynomials are given by 



= np{l — p){N + 1 — n), hn = p{N — n) + n{l — p). 



:i.6) 



Observe that = = a^jj^i- 

In this paper, we are concerned with a semi-classical generalization of Krawtchouk 
polynomials with the weight function defined by 



W[X] 



N 



X j [1 — a) 



xe{0,...,iV}, 



;i-7) 



where a < 1 and c > are two real parameters. Note that we can recover the Krawtchouk 
weight from (11.71) by letting c and —a tend to infinity and — f — ^ jz^- Our first result 
gives a discrete system satisfied by the recurrence coefficients. 



Theorem 1.1. Let an and 6„ he the recurrence coefficients in (II. 2p for the weight (II. 7p . 
We then have that 



hn + N + l + c — n — a 



N 



:i.8) 
:i.9) 



satisfy the following discrete system 



{Xn ~\~ yn){Xn+l + y n J 



{Xn + yn){Xn + 1/n-l) 

with initial conditions 

xo = 0, I/O = 



yn{N + 1 + Nyn){N +l-a + Ny^) 
cN 

Xn{-N - 1 + Nxn){a -N -1 + NXn) 



:i.io) 



N{Nxn - n) 
N + l + c-a c M(-A^ + l,2-a,-c) 



I- a M{-N, 1 - a, -c) 
where M{a,b,z) is the confluent hypergeometric function iFi{a;b;z) defined by 



[1.11) 



iFi(a; b;z) = ^ 



s=0 



{b)ss\ 



(cf. H Chapter 13]). 
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Recall that a truncated confluent hypergeometric function is a Laguerre polynomial, 



in particular M{—N, 1 — a, — c) 



-{-"), 



(l-a)jv-^Af 

value I/O is in terms of a ratio of Laguerre polynomials. 

The system (ll.lOp can be obtained by a limiting procedure from a-dPiv 
by 



— c); see |8j Eq. 13.6.19]. Hence the initial 

given 



+ Yn){Xn+l + Yn 



(Yn - A)iY^ - B){Yn - C){Yn - D) 

(y„ + r-z„)(F„-r-zo 



[1.12) 



[Xn + Yn) {Xn + Y^ 



n—l I 



{Xn + A){Xn + B){Xn + C) (X„ + D) 
{Xn + A — Zn+l/2){Xn — A — Zn+l/2) 

with A + B + C + b = Indeed, by taking 

Xn = Xn- l/e, Yn = yn+ l/s, 

A = l/€, B = -3/£+ {2N + 2- a)/N, 
C = l/e - {N +1)/N, D = l/e-{N + l-a)/N, 
Zn = {2n-l)/{2N) + l/e, T^=4c/{Ne), A = 2/e, 

and letting e tend to zero, we get fll.lOp . 

It is worthwhile to point out that if we take a = and c = jz^, the weight function 
(II. 7p reduces to the one considered in [1], where the author derived a discrete system for 
the recurrence coefficients. We also observe that the system in [3T|, which was also used in 
[To] for the generalized Meixner weight, can be obtained from fll.lOp with an appropriate 
choice of the parameters and scaling of x„ and yn- 

Since there are two parameters in the weight (11.71) . the recurrence coefficients are 
dependent on these parameters as well. Our next theorem shows that, when viewed as a 
function of c, the recurrence coefficients a„(c) and &„(c) are related to the fifth Painleve 
equation. 

Theorem 1.2. With x„(c) and yn{c) defined in (11. 8p and (II. 9p . we have 

^ ^ _ {cy' + n-N)^-2{n-N){n-N + cy')y - A^y^ + A2y' - Asy^ 



yn{c) 



AcN{y - l)?/2 

N-n + y{l + n + c- a + {a-N - l)y) - cy' 
2N{y - l)y ' 



:i.l4) 



where 

Ai = {1+N-a)^, A2 = 2{l+N-a){l+c+N-a), A3 = {l+c+n-a){l+c-n+2N-a) 
and y = y{c) is the solution of the fifth Painleve equation Py 



y 



1 

2l/ y 



1 



{y'Y - - + 



{y-^f 



B 



+ - + — + 

y 



Cy , Dy{y + l] 



y 



with the parameters given by 



B 



{n - Nf 



C = -{n + a), D = --. 



:i.i5) 



(1.16) 
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Theorem 11.21 gives the relationship between x„, ?/„ (and hence a^, 6„) and solutions 
of Py explicitly. 

The Painleve equations possess the so-called Painleve property: the solutions have 
no movable branch points. They were discovered by Painleve and his colleagues at the 
beginning of the twentieth century while classifying all second-order ordinary differential 
equations of the form 

w" = n{z,w,w'), (1.17) 

where ' = d/dz, the function TZ is rational in w and w', meromorphic in z, which pos- 
sess the Painleve property. It turns out that, up to Mobius transformations, only fifty 
equations of the form fll.l7p have the Painleve property [13], EHl [30]. Forty-four of these 
equations can either be linearized, be transformed to a Riccati equation or be solved in 
terms of elliptic functions. The six remaining equations are now known as the Painleve 
equations, which are often referred to as nonlinear special functions [6J and have numerous 
applications in mathematics and mathematical physics. 

The Painleve equations cannot be solved in terms of elementary functions or known 
classical special functions in general. For certain combinations of parameters, however, 
Pii-Pvi have solutions expressed in terms of special functions. For Py, the choice of 
parameters f ll.l6p is exactly when Py admits classical solutions expressible in terms of 
confiuent hypergeometric functions (equivalently, Kummer functions or Whittaker func- 
tions), see [HI § 32.10 (v)]. In fact, it is the case when confiuent hypergeometric functions 
have the associated Laguerre polynomials as special cases. Hence, the condition f ll.l6p 
on the parameters in Py is actually that when Py has rational solutions (see [H § 32.8 
(v)]). This is also consistent with the initial condition (11. lip for yQ. For more information 
about classical and rational solutions about Py and the associated r-functions, we refer 
to[5l[6l[T2l[l6l[I7l[ail2a[25ll271[28l[33]. 

The rest of this paper is devoted to the proofs of Theorems 11.11 and 11.21 They are given 
in Sections [2] and [3l respectively. In Section H] we study recurrence coefficients numerically. 
We conclude this paper with a discussion in Section [S] 



2 Proof of Theorem 11.1 



The proof of Theorem 11.11 relies on the ladder operators for discrete orthogonal polyno- 
mials, which were studied in [20]. We start with a brief description of this aspect. 



2.1 Ladder operators for discrete orthogonal polynomials 

Given a weight function w, we define a potential 

= _ V-(x) ^ wix-l) -wix) ^ 
w{x) w[x) 

which is a discrete analogue of the external field generated by w. The action of the forward 
difference operator 

A/(x) = /(x + l)-/(x) 
on orthogonal polynomials p„ is given by 

Ap„(x) = An{x)pn-l{x) - Bn{x)Pn{x). (2.2) 
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In case of a weight w on the lattice Nq with w{—l) = 0, the coefficients and Bn{x) 

in (12 ■2p are given by 

. / N /,N /, .u(x + l)—u(k) ,,, 

An{x) = a^Yl Pn{k)pn{k - I) \ w{k), (2.3) 

and 



Pn{k)Pn~l{k - I) ^ ' 7^w{k). (2.4) 

X + 1 — k 

feeNo 

In the case of a weight w supported on a finite lattice {0, 1, ... , A^} with boundary con- 
ditions w{—\) = and w{N + 1) = 0, the coefficients An{x) and Bn{x) in (12.21) are given 
by 

>l„(x) = + ^''"'^' .■(A^) + a„ f:p„(t)p„(* - l) ''(^ + \)-"W ^,(t). (2.5) 

fc=0 

and 

p„(iV+l)p„,i(iV) ^ 

-DnlXI Oj^ , , J- II. \- - /J- II— ± \- - -/ 

N — X ^-^ X + 

fc=0 



(iV) + a.Vp„(fc)p„_i(fc-l) "^'^ + \^ ^^^^ ^(fc). (2.6) 



In both cases, when the lattices are finite or infinite, the following compatibility relations 
between the functions An and Bn hold: 

B„{x) + Bn+i{x) = ^— ^A„(x) -u{x + l) + Y. (2-7) 



2 ^n-l(2^) 



a„+iy4„+i(x) - a„ = (x - 6„)5„+i(x) - (x + 1 - 6„)5„(x) + 1. (2.8) 

2.2 Proof of Theorem [TTI 



To prove Theorem II. we follow the theme in , which made use of the compatibility 
relations between the ladder operators. 

It is readily seen that, with w defined in (II. 7p . w{—l) = and w{N + 1) = 0. By 
(12. ip . we have 

uix) = -1 + l^l-/"^ (2.9) 
c[N + 1 — x) 

It then follows from ([23]) and ([22]) that 

Anix) (1„R„-— -|- CtnTj 



IS — X IS — X 



where 

N 



Rn = - ^^Pn{k)Pn{k - l)w{k), 



C 
k=0 
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and 

, X 1 

IS — X IS — X 

with to = and 

iV 



c 

k=0 



c A/ + 1 — /c 

fc=0 

Using the orthonormahty of the polynomials one immediately gets that 

1 

c' 

since Pn{k — 1) = p„(A;)+lower order terms. Thus, we obtain 

X 1 

■^nix) = d-n , ,r T "I" OinTi 



R = - r = n 



c{N-x) " "iV-x' 
with to = 0. 

Next we use the compatibility relations (12. 7p and f l2.8p . The first compatibility relation 
(12. 7p gives rise to the following two equations, after collecting coefficients of equal powers 
of x: 

n 

C{tn + tn+l + hnTn) = -I + a + cN + C^Tj , (2.10) 

i=o 

cTn-hn-l + a- c + n = {). (2.11) 
From the second compatibility relation we get 

al^^Tn+i - alTn-i = -6„t„+i + (6„ - l)tn + N, (2.12) 

a^+i - = c(t„,+i - in - 1). (2.13) 
From (12.131) we find, using telescopic summation, that 

al = c{tn-n). (2.14) 

Equation (12. lip gives 

hn = — c + n + a + cT„. (2-15) 

After multiplying (12.120 by T„ we can replace —hnTn using (I2.10p . which allows us to take 
a telescopic sum. As a result we get 

n— 1 71—1 

c\tn - n)T^Tn-i = ctl + cNY^ T, - cK -tn{cN-l + a). (2.16) 

j=0 i=0 



One can also replace and af^^-^^ in equation (12.121) using (I2.14p and rewrite the obtained 
equation by collecting the coefficients of and t„. Taking a telescopic sum and using 
f l?T5|) gives 

n-l 

c{tn-n){T^ + T^_,)+tr,{a + n-c-2) =nN -cY,Tj. (2.17) 

Multiplying (I2.17P by T^, using (I2.16P on the left-hand side and (I2.10p on the right-hand 
side and eliminating the sum using fl2.10p we finally obtain 

(N + cT„)(l -a-cN + ctn+i + cTn{a - c - 2 + cT„)) 

+ Ctn{N -tn+l+cTr,) =0. (2.18) 

Replacing the sum in equation (I2.17P with the aid of (I2.16P we get 

tn{l - 2n + {2n - 2 + a)iN + 1) + c{n + N)Tn + c{n + N + cT„)T„_i) 

= n{N + cT„_i)(Ar + cTn) + tl{n - 2 + a + c(T„ + T„_i)). (2.19) 

Setting 

^n = |, l/„ = -^-l (2.20) 

in the last two equations, we obtain the system fll.lOp . 

Finally, we note that the initial conditions for the recurrence coefficients, in terms of 
the moments /ii and /zq, are given by Oq = and 

/ii _ cN M{-N + 1,2- a, -c) 



/io (1 - a) M{-N, 1 - a, -c) 



where M(a, b, z) is the confluent hypergeometric function i-Fi(a; 6; z). This, together with 
( II. 8p and (II. 9p . implies initial conditions (II. lip for Xn and 
This completes the proof of Theorem 11.11 



3 Proof of Theorem [13] 

We use the method proposed in ^ [10] and system (ll.lOp to prove Theorem 11.21 We 
repeat the main steps to be self-contained. 

To prove the theorem, we derive the differential equation for In [21 [ID! we have 
used the Toda system 

K)' - 1 K) = ^(i. - 

II _ /'2 _ 2\ 

which also holds in the present case. Solving the ffist equation (ll.lOp for and the 
second equation for yn-i and substituting into the Toda system (13. ip (where we have 
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replaced and 6„ by their expressions in terms of x„ and ?/„ from Theorem II .ip . we get 
two equations 

, _ NXn-n f Xn{-N -1+ NXn){0i- N -1 + NXr. 



C \ N{NXn - n){Xn + Vn) 

and 

, _ y^i^N + 1 + Ny^){N + 1 - a + Ny,,) 

where the differentiation is with respect to c. By differentiating equation fl3.2p and sub- 
stituting the expression for x'^ we obtain an equation for y" as a function of y'^, yn, Xn- 

Vn = [NXn - n){Xn + yn) + TTT ^ ^ + Vn 

C cN[Xn + Vn) 

+ ^ V {{N + 1 + Nyn){N + 1 - a + Ny^) + iVy„(2iV + 2 - a + 2iV2/0) 

,(iV + 1 + Nyn){N + 1 - a + iVy„) 



iVc2(x„ + 7/„)2 

a;„(-iV - 1 + Nxn){a -N -1 + Nxn) 



X \Xn + yn + Cy'^-{NXn-n){Xn + yn) + 



N{Xn + yn) 



Ehminating x„ between this equation and (13.20 gives a nonhnear second order second 
degree equation for y^. 

G{ylyly^,c)=Q. (3.3) 

We have used Mathematical^] to compute this long expression. It can be checked by direct 
computations, that applying a transformation 

, , viz) + 2a - 4A^ + 2n - 3 ^ 
^"(^) = ^ ' ^ = ^ 

this equation simplifies considerably and can be written as 

{v" - Qv"^ - aiv - I3if = {v/z - 2z)^{v'^ - Av'^ - aiv^ - 2(3iv - 71), 
where v = v{z) and 

ai = 4(2a - AN + 6n - 1), 
Pi = 2{2n + l)(6n - 8A^ - 5) + 8{2n + l)a - 8a^ 
71 = 4(2a - 4A^ + 2n - 3)(4n2 + 4n + 1 - W). 

The last equation appears in [7j (equation (A. 8)) and is known to be related to the fifth 
Painleve equation. This is similar to the case of generalized Meixner polynomials (see 

0)- 

By taking y„(c) in a form as shown in (I1.14p . we get the fifth Painleve equation Py 
(ll.lSp with parameters (11.160 . 
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A similar approach can be used to obtain a (second order second degree) differential 
equation for x„. To get the expression fll.lSp for a;„(c) in terms of the solutions of the fifth 
Painleve equation we use f l3.2p . substitute (11 .Mp and find x„(c) by solving a quadratic 
equation. One can check with Mathematica that one of the roots of this equation, i.e. 
(ll.lSp . indeed gives the statement of Theorem 11.21 This completes the proof of Theorem 

Note that we can derive a nonlinear discrete second order equation for yn{c)- From 
the first equation of (ll.lOp with n and the second equation with n + 1 we eliminate Xn+i 
by computing the resultant. The obtained equation and the second equation of (ll.lOp 
with n then allow us to eliminate x„. As a result, we obtain a nonlinear discrete equation 
for yn{c) which we denote by 

FliVn-l, 2/n, Vn+l, c) = 0. (3.4) 

The equation was again obtained by using Mathematica but it is too long and too com- 
plicated to include here explicitly. Similarly, a nonlinear discrete equation for a;„(c) can 
be obtained: 

F2{Xn-l,Xn,Xn+l,c) = 0. 

One can, similarly to f?, Sect. 3], check that the functions i/n-i, Vn and Un+i in (13. 4 p are 
connected by using the Backlund transformation of the fifth Painleve equation [HI § 32.7 
(v)]. In particular, one can express all of them using only y and y' (the solution of Py 
with parameters (I1.16P ). When these expressions are substituted into (13. 4p . it becomes 
identically zero. Moreover, equation (13. 4p can essentially (up to a factor depending on 
yn) be obtained by eliminating y between 

^ -(1 + + Nyyr,){N -n + {l + N-a + NyM 
N{N -n + {l + n + NyM 

and 

^ {l + c + N-a + Ny,,-{l + N~a + Ny^P 

N{y-l){nc + N{l + c + N -a + Nyn-{l + N -a + Nyn)y)yn) ' 

where P = (A^ + l)c + N{c -n + 2 + 2N-a + Ny^ + {n-2N~2 + a- Nyn)y)yn- 




10 20 30 40 50 60 70 80 10 20 30 40 50 60 70 



Figure 1: Recurrence coefficients (a^ on the left, 6„ on the right) for generalized 
Krawtchouk polynomials (A^ = 80, a = — 1, c = 2) 
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10 20 30 40 50 60 70 }^0 H) 20 30 4U 50 60 70 80 



(a) 



(b) 



Figure 2: Recurrence coefficients (a^ on the left, 6„ on the right) for generahzed 
Krawtchouk polynomials (A^ = 80, a = 0.8, c = 2). 



4 Computing the recurrence coefficients 



1 20 30 40 50 60 70 SO 10 20 30 40 50 60 70 80 



(a) 



(b) 



Figure 3: Recurrence coefficients (a^ on the left, 6„ on the right) for generalized 
Krawtchouk polynomials (iV = 80, a = —1, c = 30). 



The system (11.101) . together with the definitions (11. 8p and (II. 9p . can be used to compute 
the recurrence coefficients starting with the initial values Xq and yo given in (II. lip . Figures 
1-3 show these coefficients for = 80 and for some different choices of the parameters 
c > and a < 1. Recall that the recurrence coefficients (II. 6p for the Krawtchouk 
polynomials are such that is quadratic in n, with Oq = = a%_^_l, and that fe„ is 
of degree 1 in n. Also recall that when c and —a tend to infinity and — ^ — ^ 
recover the Krawtchouk weight. Figure H] shows this observation. For the generalized 
Krawtchouk polynomials the are ffist increasing until they reach their maximum, and 
then decreasing, with Og = = a^+i. The 6„ are no longer linear but change from concave 
to convex. The point where reaches its maximum and the inflection point for the 6„ 
depend on the parameter c, these two points do not coincide. Moreover, for a values 
close to 1, the 6„ change once more to concave, as can be seen (albeit not very clearly) in 
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Figure 4: Recurrence coefficients (a^ on the left, 6„ on the right) for generahzed 
Krawtchouk polynomials (A^ = 80, a = —2000, c = 1000, circles) and for Krawtchouk 
polynomials (A^ = 80, p = 1/3, crosses) 



Figure |2J 

To get these values (they were obtained using Maple), one has to work with quite high 
accuracy, for it turns out that even a very small perturbation (10^^°° for = 80) in the 
initial value quickly leads to very bad results. An explanation for this behaviour is the 
following: the system f ll.lOp allows to compute the Xn, from which the coefficients can 
be obtained by fll.8p . There is, however, no a priori reason why these should be positive 
and a^^^ = 0. This is not a new observation: in [1] (Section 6.5, about discrete g-Hermite 
I polynomials) it was conjectured that the initial value needed to obtain the recurrence 
coefficients for these polynomials, is the only real value which leads to positive values for 

for every n > 0. Since in this generalized Krawtchouk case there is only a finite number 
of coefficients {of, . . . , a^) which need to be positive, this conjecture obviously does not 
hold here, but we conjecture that there is only one initial value yo for which aj^_^i = 
and the positive. An alternative, and numerically more stable way, would 

be to use a fixed point algorithm for which the fixed point satisfies fll.lOp . 

Following a referee's remark, we also used the Stieltjes approach to calculate the 
recurrence coefficients, see e.g. [HI §2.2.3.1]. In this method, one uses (11. 3p to calculate 
a„ and 6„ from pn and Pn-i- Then (II. 2p can be used to obtain an+iPn+i{x) , which after 

normalization gives Pn+i- The initial values needed are p_i = and po = , J- =. 

VSfc=o«'('=) 

Numerical calculations show that this method needs less precision in the initial conditions 
but requires more computations. 



5 Conclusion and future directions 

In this paper we have been dealing with the recurrence coefficients of the discrete orthog- 
onal polynomials, namely the generalized Krawtchouk polynomials, and have shown that 
they are related to the classical solutions (in terms of confluent hypergeometric functions) 
of the fifth Painleve equation. A possible further direction is to consider more factors in 
(12. 9 p and find out whether the sixth Painleve equation or some equation from the Painleve 
hierarchy is related to the recurrence coefficients of such a weight. A reason that leads to 
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this conjecture follows from the observation that the initial conditions for such a weight 
correspond to special solutions of the sixth Painleve equation. 
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